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Abstract 

The extrapolation of small-cluster exact-diagonalization calculations is used 
to examine the possibility of electronic ferroelectricity in the one dimensional 
spinless and spin-one-half Falicov-Kimball model (FKM). It is found that nei- 
ther spinless nor spin-one-half version of the FKM does not allow for a ferro- 
electric ground state with a spontaneous polarization, i.e. there is no nonvan- 
ishing < d + f >-expectation value for vanishing hybridization V. 
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1 Introduction 



In the last few years the FKM |1J has been extensively studied in connection with 
the exciting idea of electronic ferroelectricity 0, ^ |J. It is generally supposed that 
the ferroelectricity in mixed-valent compounds is of purely electronic origin, i.e. it 
results from an electronic phase transition, in contrast to the conventional displacive 
ferroelectricity due to a lattice distortion. Since the FKM is probably the simplest 
model of electronic phase transitions in rare-earth and transition-metal compounds 
it was natural to test the idea of electronic ferroelectricity just on this model. 

The FKM is based on the coexistence of two different types of electronic states in 
a given material: localized, highly correlated ionic-like states and extended, uncorre- 
lated, Bloch-like states. It is accepted that insulator-metal transitions result from a 
change in the occupation numbers of these electronic states, which remain themselves 
basically unchanged in their character. Taking into account only the intra-atomic 
Coulomb interaction between the two types of states, the Hamiltonian of the spinless 
FKM with hybridization can be written as the sum of four terms: 

H = E Ujdtdj + U df £ ftfidtdi + E f Y, ft fi + V £ dffi + h.c, (1) 

ij i i i 

where ft, fi are the creation and annihilation operators for an electron in the local- 
ized state at lattice site i with binding energy Ef and df , di are the creation and 
annihilation operators of the itinerant spinless electrons in the <i-band Wannier state 
at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical 
hopping of the itinerant d electrons between sites i and j. These intersite hopping 
transitions are described by the matrix elements tij, which are —t if i and j are the 
nearest neighbors and zero otherwise (in the following all parameters are measured 
in units of t). The second term represents the on-site Coulomb interaction between 
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the c?-band electrons with density rid = jJ2i dfdi and the localized / electrons with 
density n/ = ^ J2i ft fit where L is the number of lattice sites. The third term stands 
for the localized / electrons whose sharp energy level is Ef. The last term represents 
the hybridization between the itinerant and localized states. 

In spite of fact that many of ground-state properties of the spinless FKM (the 
nature of the ground state 0, the picture of valence and metal-insulator transi- 
tions m HI, etc.) are well understood at present, the problem of electronic ferroelec- 
tricity remains still an open question. Very recently Portengen at al. @, || studied the 
FKM with a k-dependent hybridization in Hartree-Fock approximation and found, 
in particular, that a no n- vanishing excitonic < d + f >-expectation value exists even 
in the limit of vanishing hybridization V — > 0. As an applied (optical) electrical field 
provides for excitations between d- and f-states and thus for a polarization expecta- 
tion value Pdf =< df fi >, the finding of a spontaneous Pdf (without hybridization 
or electric field) has been interpreted as evidence for electronic ferroelectricity. How- 
ever, analytical calculations within well controlled approximation (for Udf small) 
performed by Czycholl @] in infinite dimensions do not confirm this conclusion. In 
contrast to results obtained by Portengen et al. ||, he found that the symmetric 
(Ef = 0,nf = rid = 0.5) FKM does not allow for a ferroelectric ground state with a 
spontaneous polarization, i.e. there is no nonvanishing < d + f >-expectation value 
in the limit of vanishing hybridization. 

In order to shed some light on this controversy we have decide to study the 
problem of electronic ferroelectricity in the FKM using the small-cluster exact- 
diagonalization method that was so successful in describing of valence and metal- 
insulator transitions in this model || |6| . It should be noted that for given L the full 
Hilbert space of the spinless FKM consists of 4 L quantum states, thereby strongly 
limiting numerical computations. Although the number of states can be reduced 
considerably by the use of symmetries of H, there is still a limit (L ~ 10) on the size 
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of clusters that can be studied using small-cluster exact-diagonalization calculations. 
However, we will show later that due to small sensitivity of the FKM on L (for a 
wide range of parameters), already such small clusters can describe satisfactory the 
behavior of the < d + f >-expectation value that lies in the center of our interest. 

To compare our numerical results with Czycholl's ones M obtained for small 
Coulomb interactions we have started our investigation in the weak-coupling limit 
and half-filled band case (n/ = rid = 0.5). The weak-coupling numerical results for 
Pdf obtained using the modified Lanczos method f|, || are displayed in Fig. 1. To 
reveal the finite-size effects on numerical calculations have been performed for 
three finite clusters of L = 6, 8 and 10 sites. It is seen that there are nonzero finite- 
size effects on the < d + f >-expectation value in the weak-coupling limit, however 
they do not change qualitatively the behavior of Pdf in the limit V — > that is crucial 
for the verification of spontaneous polarization. In all cases the < d + f >-expectation 
value vanishes in the limit V — > 0, so there is no spontaneous polarization in the 
spinless FKM. Thus, in accordance with Czycholl's weak-coupling results we can 
conclude that the spinless FKM does not allow for a ferroelectric ground state with 
a spontaneous polarization at least for Udf small. 

Unlike the method used by Czycholl that is restricted to small interactions we can 
proceed in the numerical study of the FKM at arbitrary Udf- The strong-coupling 
numerical results for Pdf are displayed in Fig. 1 for Udf = 4 and Udf = 10. Obvi- 
ously the one-dimensional FKM does not exhibit a ferroelectric ground state with 
a spontaneous polarization in the strong-coupling limit. For both values of Udf the 
< d + f >-expectation value vanishes for V — > 0, and it is demonstrated that this 
result is independent of L. Thus, the strong-coupling results can be satisfactory 
extended to large systems and they should be considered as definite. 

Of course, the absence of ferroelectric ground state in the spinless FKM does 
not exclude that some other models could exhibit such a ground state. The spinless 
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FKM is not too realistic model of a rare-earth compound, because any real Fermi 
system has at least a spin degeneracy. Therefore, it is natural to ask if the spin- 
one-half FKM would not allow for the ferroelectric ground state with a spontaneous 
polarization. 

Numerically the problem can be easily solved since the numerical method used 
for the spinless FKM can be straightforwardly generalized also for the spin-one- 
half FKM. Unfortunately, including spins will result in further reduction of the size 
of clusters that can be analyzed using the exact-diagonalization method. In order 
to compensate partially for the small size of clusters we next examine the model 
only for strong d-f interactions that (as was shown for the spinless FKM) minimize 
considerably the finite-size effects. 

The Hamiltonian of the spin-one-half FKM can be obtained directly from the 
spinless model by including the spins for both d and / electrons and by adding the 
on-site Coulomb interaction Uff that acts between two / electrons of opposite spins 
(the last term): 



H = ^t ij d&d j „ + Uf^f&f u ,d&d u , + E f ^f£f u , + V^d£f Ur + h.c. 

icrcr' icr i(T 

fiaficrfi-afi—cr- (2) 



i j a 



The ground-state properties of this model, for V = 0, have been investigated in 



our preceding paper | 10| . We have found that numerical results depend strongly on 
/-/ interaction strength Uff, but they are relatively insensitive to d-f interaction 
strength U^f (for U^f > 2). Therefore to represent the typical behavior of the model 
at nonzero V, and to minimize finite-size effects we choose in the next study the 
value Udf = 3 that is sufficiently large to stabilize the system. Another advantage 
of this selection is that the ground-state phase diagram of the spin-one-half FKM 
without hybridization is well understood for large values of Udf- Particularly, in 



the strong-interaction limit Uff > 4/7T the ground state is insulating for Ef < —A/tt 
and metallic for Ef > —A/tt. At Ef = —4/tt the model exhibits a discontinuous 
insulator-metal transition that is accompanied by a integer-valence transition from 
rif = 1 (n<i — 0) to rif = (rid = 1). The same behavior exhibits the model also in 
the opposite limit Uff < 4/tt: the ground state is insulating for Ef < E c (Uff) and 
metallic for Ef > E c {Uff). However, a discontinuous insulator-metal transition that 
takes place at Ef = E c (Uff) realizes now between an integer-valence state rif = 1 and 
an inhomogeneous intermediate- valence state rif ^ 0. These results show that there 
are only three physically different ground states in the spin-one-half FKM without 
hybridization, and namely, an insulating integer- valence ground state with rif = 1, 
a metallic integer- valence ground state with rif = 0, and a metallic intermediate- 
valence ground state with < rif < 1. Here we examine whether these ground states 
are stable against a small, finite hybridization or whether some new ground states 
are obtained if one starts from a finite hybridization and studies the V — > limit of 
the model. Again the special attention is devoted to the question if the model can 
exhibit a ferroelectric ground state with a spontaneous, nonvanishing polarization 

Pdf — < die fia >. 

The strong-coupling [Uff = 10) numerical results for Pdf obtained on finite clus- 
ters of 4 and 6 sites are displayed in Fig. 2. It is seen that for both Ef < — 4/-7T (for 
V = an insulating integer- valence state) and Ef > — 4/-7T (for V = a metallic 
integer-valence state) the < d + f >-expectation value vanishes in the limit V 
indicating that there is no spontaneous polarization in the spin-one-half FKM. It 



should be noted that this result is expected since approximate solutions |TT| lead 
also to Pdf = for the integer- valence states with rif = 1 and rif = 0. A less trivial 
situation is expected in the intermediate- valence state with < rif < 1. As was 
discussed above this state exists in the spin-one-half FKM without hybridization for 
Uff < 4/tt and Ef > E c (Uff). Particularly, for L = 4, U^f = 3 and Uff = 0.4 we have 
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found that n f = 1 for E f < -1.35, n f = 0.5 for -1.35 < E f < -0.65, and n f = 
for > —0.65. To examine the stability of intermediate-valence state against a 
small, finite hybridization we chose the value Ef — — 1 for numerical calculations. 
The results obtained for P d f are shown in Fig. 3. Again the < d + f >-expectation 
value vanishes in the limit V — >• indicating the absence of spontaneous polarization. 
We have verified this important result also for L = 8 (see inset in Fig. 3). Unfortu- 
nately, due to the memory limitations we were not able to continue with numerical 
computations on larger lattices (the memory requirement of the Lanczos method for 
clusters larger than L = 8 is beyond the reach of present day computers) to exclude 
definite the existence of spontaneous polarization in the spin-one-half FKM. How- 
ever, in accordance with results obtained for the spinless FKM we do not expect that 
the behavior of Pjf for V — >• could be qualitatively changed on larger lattices. 

In summary, we have used the extrapolation of small-cluster exact-diagonalization 
calculations to study the possibility of electronic ferroelectricity in the one dimen- 
sional spinless and spin-one-half FKM. It was found that neither spinless nor spin- 
one-half version of the FKM does not allow for a ferroelectric ground state with a 
spontaneous polarization, i.e. there is no nonvanishing < d + f >-expectation value 
for vanishing hybridization V. 

This work was supported by the Slovak Grant Agency for Science under grant 
No. 4177/97. 
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Figure Caption 



Fig. I. Hybridization dependence of the d- /-polarization P d f =< dj f { > in the 
spinless FKM calculated for three different values of Uff and L. The symmetric case 
E f = 0. 

Fig. 2. Hybridization dependence of the d- /-polarization P d f =< df a f icr > in the 
spin-one-half FKM calculated for Ef = —2 (for V = an insulating integer- valence 
state) and Ef = (for V = a metallic integer- valence state). Inset: Hybridization 
dependence of P d f in the limit of vanishing V for L = 4 and L = 6. 

Fig. 3. Hybridization dependence of the d- /-polarization P d f =< df a f ia > in 
the spin-one-half FKM calculated for Ef = —1 (for V = an intermediate-valence 
state). Inset: Hybridization dependence of P d f in the limit of vanishing V for L = 4 
and L — 8. 
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